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Singular limits for 4-dimensional semilinear elliptic problems 
cu with exponential nonlinearity 
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.^^ • Abstract 

C) ■ Using some nonlinear domain decomposition method, we prove the existence of singular 

limits for solution of semilinear elliptic problems with exponential nonlinearity. 

1 Introduction and statement of the results 



> 

in 



In the last decade important work has been devoted to the understanding of singularly perturbed 
problems, mostly in a variational framework. In general, a Liapunov-Schmidt type reduction 
argument is used to reduce the search of solutions of singularly perturbed partial differential 
equations to the search of critical points of some function which is defined over some finite 
C^) ' dimensional domain. 

One of the purpose of the present paper is to present a rather efficient method to solve such 
^sO . singularly perturbed problems. This method has already been used successfully in geometric 

context (constant mean curvature surfaces, constant scalar curvature metrics, extremal Kahler 
metrics, manifolds with special holonomy, . . . ) but has never appeared in the context of partial 
differential equations. We felt that, given the interest in singular perturbation problems, it was 
worth illustrating this on the following problem : 
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Let il C M 4 be a regular bounded open domain in R 4 . We are interested in positive solutions 

of 
K/ i f A 2 u = pV in fi 

(1) 
[ u = Au = on dfl, 

when the parameter p tends to 0. Obviously, the application of the implicit function theorem 
yields the existence of a smooth one parameter family of solutions (u p ) p which converges uni- 
formly to as p tends to 0. This branch of solutions is usually referred to as the branch of 
minimal solutions and there is by now quite an important literature which is concerned with the 
understanding this particular branch of solutions ^3] . 

The question we would like to study is concerned with the existence of other branches of 



solutions as p tends to 0. To describe our result, let us denote by G(x, •) the solution of 

f A 2 G(a;,-) = 64 ir 2 6 X in fl 

(2) 
1 G(x,-) = AG(av) = on 5JL 

It is easy to check that the function 

R(x,y):=G{x,y) + 8log\x-y\ (3) 

is a smooth function. 
We define 

m 

W(x\...,x m ) :=^i?(^',a; J ) + ^G(x J ,^). (4) 

Our main result reads : 

Theorem 1.1 Assume that (x , . . . , x™ 1 ) is a nondegenerate critical point ofW, then there exist 
po > anrf (u p )pg(o iPo ) a one parameter family of solutions of £|J), smc/j i/iat 



limu p = ^G(a; J ,-) 



mClO-R...,,'"}). 

This result is in agreement with the result of Lin and Wei [7| (see also ^Sl) where sequence of 
solutions of Q) which blow up are studied. Indeed, in this paper, the authors show that blow up 
points can only occur at critical points of W. 

Our result reduces the study of nontrivial branches of solutions of Q to the search for critical 
points of the function W defined in |0J. Observe that the assumption on the nondegeneracy of 
the critical point is a rather mild assumption since it is certainly fulfilled for generic choice of the 
open domain Q. 

Semilinear equations involving fourth order elliptic operator and exponential nonlincarity 
appear naturally in conformal geometry and in particular in the prescription of the so called 
Q-curvature on 4-dimensional Riemannian manifolds |3], 0] 

Q 9 = ^(-A 9 S s + S 3 2 _ 3 |Ric 9 | 2 ) 

where Ric g denotes the Ricci tensor and S g is the scalar curvature of the metric g. Recall that 
the Q-curvature changes under a conformal change of metric 

g w = e 2w g, 



according to 
where 



P g w + 2Q g = 2Q gvj e 4w (5) 



P g :=A a g + 6(jS g I-2Bic g ),l «>i 



is the Panietz operator, which is an elliptic 4-th order partial differential operator 4 and which 
transforms according to 

e Aw P e ^ g = P g , (7) 

under a conformal change of metric g w := e 2w g. In the special case where the manifold is the 
Euclidean space, the Panietz operator is simply given by 

P -A 2 

in which case (JHJ reduces to 

A 2 w = Qe 4w 

the solutions of which give rise to conformal metric g w = e 2 w g e ud whose Q-curvature is given 
by Q. There is by now an extensive literature about this problem and we refer to 0] and [H] for 
references and recent developments. 

When n = 2, the analogue of the Q-curvature reduces is the Gauss curvature and the corre- 
sponding problem has been studied for a long time. More relevant to the present paper is the 
study of nontrivial branches of solutions of 

-Au = p 2 e u in fi 

(8) 
u = on <9f2, 

which are defined on some domain of R 2 . The study of this equation goes back 1853 when 
Liouville derived a representation formula for all solutions of (jHJ which are defined in R 2 , jTU] . 

It turn out that, beside the applications in geometry, elliptic equations with exponential non- 
linearity also arise in the modeling of many physical phenomenon such as : thermionic emission, 
isothermal gas sphere, gas combustion, gauge theory |17) . . . . 

When p tends to 0, the asymptotic behavior nontrivial branches of solutions of (JSJ is well 
understood thanks to the pioneer work of Suzuki [16| which characterizes the possible limit of 
nontrivial branches of solutions of JHJ ■ The existence of nontrivial branches of solutions was first 
proven by Weston J2HJ and then a general result has been obtained by Baraket and Pacard |2]. 
More recently these results were extended, with applications to the Chern-Simons vortex theory 
in mind, by Esposito, Grossi and Pistoia [H] and Del Pino, Kowalczych and Musso [S] to handle 
equations of the form 

-Au = p 2 Ve u 



where V is a non constant (positive) function. We give in Section 9 some results concerning the 
fourth order analogue of this equation. Let us also mention that the construction of nontrivial 
branches of solutions of semilinear equations with exponential nonlinearities has allowed Wente 
to provide counterexamples to a conjecture of Hopf 19 concerning the existence of compact 
(immersed) constant mean curvature surfaces in Euclidean space. 

We now describe the plan of the paper : In Section 2 we discuss rotationally symmetric 
solutions of J3>. In Section 3 we study the linearized operator about the radially symmetric 
solution defined in the previous section. In Section 4, we discuss the analysis of the bi-Laplacc 
operator in weighted spaces. Both section strongly use the b-operator which has been developed 
by Melrose jTJ] in the context of weighted Sobolev spaces and by Mazzeo m] in the context of 
weighted Holder spaces (see also |TS]'). 

A first nonlinear problem is studied in Section 6 where the existence of an infinite dimensional 
family of solutions of Q which are defined on a large ball and which are close to the rotationally 
symmetric solution is proven. In Section 7, we prove the existence of an infinite dimensional 
family of solutions of QJ which are defined on fl with small ball removed. Finally, in Section 8, 
we show how elements of these infinite dimensional families can be connected to produce solutions 
of Q) described in Thcorcm ll.il This last section borrow ideas from applied mathematics were 
domain decomposition methods are of common use. In Section 9 is devoted to some comments. 
In Section 10, we explain how the results of the previous analysis can be extended to handle 
equations of the form A 2 u = p 4 V e u . 

2 Rotationally symmetric solutions 

We first describe the rotationally symmetric solutions of 

A 2 u - p A e u = 0, (9) 

which will play a central role in our analysis. By the classification given by 6 , for e > 0, we 
define 

u e (x) := 4 log(l + e 2 ) - 4 log(e 2 + |a;| 2 ). 

which is clearly a solution of © when 
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(1 + e 2 ) 4 ' 



Let us notice that equation © is invariant under some dilation in the following sense : If u is 
a solution of © and if r > 0, then u(t •) + 4 log t is also a solution of 10. With this observation 
in mind, we define, for all r > 

u s , T (x) := 4 log (1 + e 2 ) + 4 log r - 4 log (s 2 + t 2 \x\ 2 ). (11) 



3 A linear fourth order elliptic operator on R 4 

We define the linear fourth order elliptic operator 

which corresponds to the linearization of © about the solution U\(= u s —i) which has been 
defined in the previous section. 

We are interested in the classification of bounded solutions of L w = in R 4 . Some solutions 
are easy to find. For example, we can define 

1-r 2 
</>oO) :=rd r ui(x) +4 = 4 — — j, 

1 + r A 

where r = |x|. Clearly h<fi — and this reflects the fact that J5J) is invariant under the group of 
dilations r — ► u(t ■) + 4 log r. We also define, for i = 1, . . . , 4 

<f>i(x) := -d Xi Ui(x) = - l 2 , 
1 + \x\ z 

which are also solutions of L <j>j — since these solutions correspond to the invariance of the 
equation under the group of translations a — > u(- + a). 

The following result classifies all bounded solutions of L w = which are defined in R 4 . 

Lemma 3.1 Any bounded solution o/Lio = defined in ]R 4 is a linear combination of 4n for 
* = 0,1,...,4. 

Proof : We consider on R 4 the Euclidean metric g euc i = dx 2 and the spherical metric 

9si = (1+tl 2 ) 2 ^ 

induced by 

II: R 4 — > S 14 

T . , ( 2x 1-I^H 

the inverse of the stereographic projection. 

According to |3] we have P g 4 = A| 4 — 2 A52 and -P gellci = A 2 . Therefore, we obtain from J7J 

^ fA 2 c4 -2Ac2) = A 2 

(l + |a;| 2 ) 2 y { s S ' 



In particular, if w : R 4 — ► R is a bounded solution of Lw = then, w : S 4 — {N} — ► R is a 
bounded solution of 

(A| 4 - 2 A S 4 - 24) w = (13) 

away from the north pole N £ S 4 . It is easy to check that the isolated singularity at the north 
pole is removable (since w is assumed to be bounded) and hence (|13[1 holds on all S 4 . 

We now perform the cigcnfunction decomposition of w in terms of the eigendata of the Lapla- 
cian on S 4 . We decompose 

w = y. w i 

e>o 

where we belongs to the £-th eigenspace of A 54, namely, we satisfies Astute — —Xewe with 

X e :=£(£ + S). 

We get from (jT3|) 

(Af + 2 X t - 24)^ = 0. 

Hence, iu^ = for all £ except eventually those for which Xg = 4. This implies that w : S 4 — ► R 
is a combination of the eigenfunctions associated to £ = 1 which are given by y>i(y) = j/i for 
i = 1, . . . , 5, where y — (yi, . . . , j/5) £ S* 4 . The sphere being parameterized by the inverse of the 
stereographic projection we may write y = U(x). Then, the functions Aipi precisely correspond 
to the functions (j>i for i = 1, . . . ,4, while the function 4^5 corresponds to the function <^o- This 
completes the proof of the result. □ 

Let B r denote the ball of radius r centered at the origin in R 4 . 

Definition 3.1 Given k £ N, a £ (0, 1) and fi £ R, we define the Holder weighted space C^ ,Q (R 4 ) 
as the space of functions w £ C /o '"(R 4 ) for which the following norm 

IHIcfT (R4) := IMIc*.»(Si) + SU P ((! + r 2 )-^ 2 \\w{r OllcJ-f^-Bwa)) > 

is finite. 

More details about these spaces and their use in nonlinear problem can be found in 15 . Roughly 
speaking, functions in C^'°(R 4 ) are bounded by a constant times (1 + r 2 ) M / 2 and have ^-th partial 
derivatives which are bounded by (1 + r 2 )^ - ^' 2 , for £ — 1, . . . , k + a. 

As a consequence of the result of Lemma 13.11 we have the : 
Proposition 3.1 Assume that \i > 1 and n $■ N, then 

L,: C 4 ^(R 4 ) — C^ 4 (R 4 ) 
w 1 — ► Ltd 

is surjective. 



Proof : The mapping properties of L M are very sensitive to the choice of the weight fi. In 
particular, it is proved in J2] and ^J (see also ^3) that L^ has closed range and is Fredhlom 
provided fi is not an indicial root of L at infinity. Recall that £ £ R is an indicial root of L at 
infinity if there exists a smooth function v on S 3 such that 

L(|x|^) = e>(|z| f - 5 ) 

at infinity. It is easy to check that the indicial roots of L at infinity are all (, £ Z. Indeed, let e 
be an eigenfunction of A53 which is associated to the eigenvalue 7 (7 + 2), where 7 € N, hence 

A53 e = — 7 (7 + 2)e. 

Then 

L (|x|« e) = (C - 7)(C - 7 - 2) (C + 2 + 7) (C + 7) M C ~ 4 e + 0(\x\^ s ). 

Therefore, we find that —7—2, —7, 7 and 7+2 are indicial roots of L at 0. Since the eigenfunctions 
of the Laplacian on the sphere constitute a Hilbert basis of L 2 (S 3 ), we have obtained all the 
indicial roots of L at infinity. 

If (i (fc Z, some duality argument (in weighted Sobolev spaces) shows that the operator L^ is 
surjective if and only if the operator £_ M is injective. And, still under this assumption 

dimKcrL M = dimCokcr L_ M . 

The result of Lemma lrt.ll preciselv states that the operator L^ is injective when jj, < — 1. Therefore, 
we conclude that L M is surjective when /j, > 1, \i $ Z. This completes the proof of the result. □ 

4 Analysis of the bi-Laplace operator in weighted spaces 

Given x 1 , . . . , x m £ fl we define X := (x 1 , . . . , x m ) and 

fl*(X) :=ti~{x\...,x m }, 

and we choose r$ > so that the balls B ro {x l ) of center x l and radius r$ are mutually disjoint 
and included in $7. For all r £ (0, ro) we define 

firPQ ■=n-{j'jL 1 B r (x j ) 

With these notations, we have the : 

Definition 4.1 Given k £ R ; a £ (0, 1) and v £ W, we introduce the Holder weighted space 
C^ ,a (ti* (X)) as the space of functions w £ C ;o '"(f2* (X)) which is endowed with the norm 

m 

IMIc*.«(fi. (x)) '■= IMIc*><*(n ro/2 (x)) + Y. SU P ( r ~ V IM** + r Olle^^-so) > 

j=l r€(0 : r /2) 

is finite. 



Again, these spaces have already been used many times in nonlinear contexts and we refer to 
|T5] for further details and references. Functions which belong to C** a (fi* (X)) are bounded by 
a constant times the distance to X to the power v and have £-th partial derivatives which are 
bounded by a constant times the distance to X to the power v — £, for i — 1, . . . , k + a. 

When k > 2, we denote by [C*' a (fi* (X))] be the subspace of functions w G C*' a (fi* (X)) 
satisfying u) = Aw = 0. 

We will use the following : 
Proposition 4.1 Assume that v < and v ^1, then 

C v :[#«(fi*(*))]o — > C°L a 4 (fi*(X)) 



is surjective. 

Proof : Again this result follows from the theory developed in ^2] an d EH ( see a ^ so USD- The 
mapping properties of £„ depend on the choice of the weight v. The operator £„ has closed range 
and is Fredhlom provided v is not an indicial root of A 2 at the points & . Recall that £ G M is 
an indicial root of L v at x J if there exists a smooth function v on S* 3 such that 

h(\x-X j \ c v)=0(\x-X j \ c - 3 ) 

at a;- 7 . As in Proposition 14. II it is easy to check that the indicial roots of C at £■? are all C, G Z. 

If ;/ ^ Z, some duality argument (in weighted Sobolev spaces) shows that the operator C v is 
surjective if and only if the operator £_„ is injective. And, still under this assumption 

dimKcr£„ = dimCokcr£_, y . 

We claim that the operator £„ is injective if v > 0. Indeed, isolated singularities of any 
solution w G C^'"(n* (X)) of A 2 u> = in Q* are removable if v > 0. Therefore, w is a bi- 
harmonic function in £1 with it) = Aw = on <9£7. This implies that w = and hence £„ is 
injective when v > as claimed. 

We then conclude that C v is surjective when v < 0, i/ ^ Z. This completes the proof of the 
result. □ 

Given y 1 , . . . , y m close enough to x 1 , . . . , x m , we set Y := (y 1 , . . . , y m ) and we define a family 
of diffeomorphisms D(— Dx.y) 

D-.n — >n 

depending smoothly on y 1 , . . . , y m by 

m 

D(s) = x + ^ xr (a; - x J ) (x J - ^'), (14) 



where Xr is a cutoff function identically equal to 1 in i? ro /2 & n d identically equal to outside 
B ro . In particular, D{y 3 ) — x° for each j. 

The equation A 2 w = f where / E C v " A {Vt(Y)) can be solved by writing w = w o D and 
/ = / o D so that w is a solution of the problem 

A 2 w+ (A 2 (woD)-{A 2 w)oD) oD- 1 = f (15) 

where this time / e C°'"(ft (X)). It should be clear that 

|| (A 2 (w o £>) - (A 2 w) o £>) o D- 1 || c o.- ((fi . (x)) < c||Y - X|| |Mb« ( fi. (JC)) (16) 

provided ||F-X|| < r /4. 

We fix ^ < 0, v $■ Z and use the result of Proposition 14.11 to choose a right inverse Gv,x 
for £„ : [C 4 > Q (iT (X))] — > C£? 4 (fi* (X)). The estimate JTBJ together with a perturbation 
argument, shows that (|15|l is solvable provided Y is close enough to X. This provides a right 
inverse Q v y which depends continuously (and in fact smoothly) on the points y , . . . , y m in the 
sense that 

/ >— G^YifoDx^oiDx.y)- 1 

depends smoothly on Y . 

5 Bi-harmonic extensions 

Given tp e C 4 > a (S 3 ) and tp E C 2 ' a (S 3 ) we define W{= W(tp, tp ; •)) to be the solution of 

A 2 W = in B 1 

W = tp on dB x (17) 

AW = ip on dB u 

where, as already mentioned, Bi denotes the unit ball in M 4 . 

We set B* = B\ — {0}. As in the previous section, we define : 

Definition 5.1 Given k £ N, a G (0, 1) and /i G K, we introduce the Holder weighted spaces 
C^ a {B1) as the space of function in C l( ^"(Bl) for which the following norm 



IMIdrfsn = su p ( r "IM»"Ollc*.«(fl a -.Bo)' 

M x r<l/2 



is finite. 



This corresponds to the space and norm already defined in the previous section when Q = B\, 
m = 1 and x 1 = 0. 

Let ei, . . . , e4 be the coordinate functions on S 3 . We prove the : 



Lemma 5.1 Assume that 



8 <p — ip) dvss = and also that / (12 if — ip) eg dv$3 = (18) 

s 3 ' Js 3 

for £ = 1, . . . , 4. TTien i/iere exists c > swc/i i/iai 

ll^fa. V> ; Ollc^BI) ^ c (IMIc«.«»(s») + IMIc«(S»))- 



Proof : There are many ways to proof this result. Here is a simple one which has the advantage 
to be quite flexible. We consider the eigenfunction decomposition of ip and ip in terms of the 
cigcnfunctions of Ag3 . 

<p = y]<Pt and i> = y~]'<Pe, (19) 

e>o e>o 

where, for each £ > 0, the functions tpg and ^ belong to the ^-th cigenspace of A53, namely 

A S 3ip e = -£{2 + £)ipg and A S 3tf> e = -£ (2 + £) ipg. 
Then the function H l can be explicitly written as 

"'-E^-i^V^ii^V"*' (20) 

Observe that, under the hypothesis, the coefficients of r° and r 1 vanish and hence, at least 
formally, the expansion of H only involves powers of r which are greater than or equal to 2. 

We claim that 

||<^||l« ^ c i IMU 2 IIVvlU- < c e \\ip\\ L 2 

where the constant eg depends polynomially on £. For example, we can write ipg = ag eg where 
ag £ R and eg is an eigenfunction of Aga which is normalized to have L 2 norm equal to 1. Then 



ipe eg dv S 3 
s 3 



< cWvWl 2 < c|M|z,° 



Next, eg solves Agseg — — £(2 + £)eg, we can use elliptic regularity theory to show that the 
L°°(S 3 ) norm of eg depends polynomially on £. The claim then follows at once. 

This immediately yields the estimate 

sup (r- 2 |JET*| + \AW\) < c(\\<p\\ L °o + ||V>|U~) 

r<l/2 

This estimate, together with the maximum principle and standard elliptic estimates yields 

supr- 2 |ir|<c(|M| L oo + ||,/;|| Loo ) 
r<l 

10 



The estimate for the derivatives of H l now follows at once from Schauder's estimates. □ 

Given 93 <E C 4 - a (S 3 ) and ip £ C 2 > a (S 3 ) we define (when it exists !) H e (= H e (ip,?P ;■)) to be 
the solution of 

A 2 H e = in K 4 - B x 

H e = Lp on dB x (21) 

AH e = V on dB t , 
which decays at infinity. 

Definition 5.2 Given k £ N, a £ (0, 1) and v £ M., we define the space C£" ,Q (R 4 — B{) as the 
space of functions w £ C ; j"(R 4 — B\) for which the following norm, 

IMIc*-(R«-.BO =SU P ^""Ik^OllcJ'-ffia-Bi)) 1 

is finite. 

We prove the : 

Lemma 5.2 Assume that 

I ipdv S 3 =0. (22) 

Js 3 

Then there exists c > smc/i t/iat 

ll-^OP.^OIIci-fOR'-Bx) ^ c (llvllc*.«(S3) + IIV'llc^CSS))- 



Proof : We use the notations of the previous Lemma. Now, the function H e can be explicitly 
written as 

H> = r- ^ + £ r" 9 ^ (w + ll A ~ E A ^ ** ^ 

e>i ^ ' £>i 

Observe that, l|22|l implies that the coefficients of r°, vanishes and hence the expansion of H e 
only involves powers of r which are lower than or equal to —1. The proof is now identical to the 
proof of Lemma r5.1l and left to the reader. □ 

Under the hypothesis of the Lemma 15.11 there is uniqueness of the bi- harmonic extension of 
the boundary data which decays at infinity. 

If F C L 2 (S 3 ) is a space of functions defined on S 3 , we define the space F to be the subspace 
of functions of F which arc L 2 (S 3 )-orthogonal to the functions 1, ei, . . . , e^. We will need the : 



11 



Lemma 5.3 The mapping 

V: C^t^xC'^S 8 ) 1 — > ^"(SfxC^fS 3 ) 1 

(v,VO i — ► {d r H l - d r H e , d r AW - d r AH e 

where H l = H l (ip, ip ; •) and H e = H e ((p, ip ; ■), is an isomorphism. 
Proof : Granted the explicit formula given in the previous two Lemmas, we have 



(24) 



We denote by W k,2 (S 3 ) Sobolev space of functions whose weak partial derivatives, up to order k 
are in L 2 (S 3 ). The norm in W k ' 2 (S 3 ) can be chosen to be 



1/2 
L 2 (S 3 ) 



\\<p\\w..* ( s*) '•= £(i + *) afc llw" 2 

\e>o 
when the function </? is decomposed over eigenspaces of A S 3 

e>o 

where A53 ipi — — £ (£ + 2) ip£. It follows at once that 

V : W k+3 ' 2 (S 3 ) X x W k+l ' 2 (S 3 ) x — * W k+2 ' 2 (S 3 ) X x W k ' 2 (S 3 ) x 

is invertible. Elliptic regularity theory then implies that the corresponding map is also invertible 
when defined between the corresponding Holder spaces. □ 

6 The first nonlinear Dirichlet problem 

For all e, r > 0, we set 



Re := T fy/i. 

Given ip G C 4 > a (S 3 ) and if) E C 2 - a (S 3 ) satisfying JTSJl, we define 

u:= Ul +H l (ip,^; (-/R £ )). 
We would like to find a solution u of 

A 2 u-24e n = (25) 

12 



which is defined in Br s and which is a perturbation of u. Writing u = u + v, this amounts to 
solve the equation 

Lv = ^ay^ WM ' A -' R ' ))+V " 1 " «). (26) 

since H l is bi-harmonic. 

We will need the following : 

Definition 6.1 Given f > 1, k € N, a € (0, 1) and /i € M, i/ie weighted space C^' a (Bf) is defined 
to be the space of functions w € C k ' a (Bf) endowed with the norm 

; llc"' a (s,) := \\ w \\c"--(B 1 ) + sup (r^ \\w(r •)IIc«'.<'(b 1 -b 1/3 )) ■ 



w 



For all ct > 1, we denote by 
the extension operator defined by 

M/)M = x(M)/(^ 

where t i — ► x(t) is a smooth nonnegative cutoff function identically equal to 1 for t > 2 and 
identically equal to for t < 1. It is easy to check that there exists a constant c = c(fi) > 0, 
independent of a > 1, such that 

||£ CT (w)|| c o, Q(R4) < c|H| c o,« (Sct) . (27) 

We fix 

/xe (1,2) 

and denote by Q^ a right inverse provided by Proposition 13. II To find a solution of (|26|) . it is 
enough to find u e C 4 ' Q (R 4 ) solution of 

v = N(e,r,<p,il>;v) (28) 

where we have defined 

N(e,r,^;v):=g,oS Re ( ^^ w (gHW, (■/*.))+« - l - „)) 

Given k > 1 (whose value will be fixed later on), we now further assume that the functions 
<p e C 4 - Q (5 3 ), -0 € C 2 - Q (S' 3 ) and the constant r > satisfy 

|log(r/r*)| < /telogl/e, ||^|| C 4,« ( <y 3 ) < «e and ||V'llc».<«(53) < « £ > ( 29 ) 

where T* > is fixed. 

We have the following technical : 
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Lemma 6.1 Given k > 0. There exist e K > 0, c K > and c K > such that, for all s £ (0,e K ) 



\\N(£,T,ip,ilj; 0)|| C 4,« (K4) < c K i 



(30) 



Moreover, 



and 



\\N(e,T,tp,ip; v 2 ) - N(e,T,cp,ip; Ui)|| c *.« (M 4) < c K e 2 \\v 2 - «i|| c *.« (K 4) 

||iV(e, r, ^ 2 , 1P2 ; v) - N(s, t, ipi,ipi ; w)|| c *.« (R 4) 

< c K e 3 (H^-^illc^-fs 3 ) + Il^2-^i||c 2 -(s3)) 
provided v = v, vi, V2 € C 4,a (IR 4 ), ip = cp, tpi, ip 2 £ C 4,Q (S' 3 ) 7 ip — ip, ipi, ip 2 € C 4,Q (S' 3 ) satisfy 



(31) 



(32) 



||W|| C 4, Q(K4) < 2C K £ 2 , ||0|| C 4,<»( S S) < KE, ||^|| C 2, Q(S 3) < K£, 

and |log(r/r*)| < /celogl/e. 

Proof : The proof of these estimates follows from the result of Lemma 15.11 together with the 
assumption on the norms of <p and i/j. Let c£ denote constants which only depend on k (provided 
e is chosen small enough). 

It follows from Lemma 15. II that 

IliTV^; ■/R E )\\cf{B RE ) < cR e 2 (IMIc*.<.(s») + U\\c^{S*)) <ci 1] e 2 
Therefore, we get 

< c^ e 2 



(1 + |- | 2 )- 4 ( e "W;7**)_i) 



Making use of Proposition 13 . II together with i|27|) we conclude that 

\\N(e,T,tp,ip; 0)|| C 4, a(K4) <c K e 2 . 
To derive the second estimate, we use the fact that 



(1 + I • | 2 )- 4 e ^(v^;-/fl e ) ( e « _ e *i _ W2 + Wl ) 



and 



(l + H 2 r 4 (e"^-/^-l) (va-t*) 
provided wi,«2 € C 4 ' a (R 4 ) satisfy ||^i|| c 4,<*, E4 * < 2c K e 2 . 



<C( 3 ) £ 2 ||« 2 - V1 || C 4, Q 



<4 4) e 2 ||U2-Ul|| C 4,a (R4) , 
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Finally, in order to derive the third estimate, we use 

(1 + | . | 2 )- 4 ( e H*(v2^2 ; -/Re) _ e H%i,-0i ; -/Re)\ £ v 



c°i Q 4 (£«J 



< cl 5) e 2 11^(^2 - <pi, fa - 4>i ; ■/R e )\\ c *,« iSitt) , 

provided v £ C 4,Q (1R 4 ) satisfies ||u|| t ; 4 ' t *(M 4 ') — 2c K e 2 . The second and third estimates again follows 
from Proposition 13.11 and (|27[) . □ 

Reducing e K > if necessary, we can assume that, 

~^e 2 <\ (33) 

for all e £ (0, e K ). Then, (|3L)|) and i|31|) in Lemma lb. II are enough to show that 

vi — ► N(e,T,<p,i>; v) 

is a contraction from 

{« € C*' a (R 4 ) : |M| c 4,„ (r4) <2c k£ 2 } 

into itself and hence has a unique fixed point v(e, r, </?,?/' ; •) in this set. This fixed point is a 
solution of (|26|l in Bn m ■ 

Reducing e K if this is necessary, it follows from (|31fl and l|32l) in Lcmma lixll tliat 
\\v(e,T,ip2,ip 2 ; •) -v(e,T,<pi,ipi ; -)ll < 2c K e 3 (\\<p 2 - ¥>i||c*.<»(S3) + 11^2 - Vi||c 2 -°(s 3 )) • (34) 
We summarize this in the : 

Proposition 6.1 Given k > 1, i/iere exist e K > and c K > (WZy depending on k) such that 
given tp € C 4 ' Q (S' 3 ) ; ^ € C 2 '"(S' 3 ) and r > satisfying fity and 

|log(r/r*)| < Acelogl/e, ||^|| C 4, Q(S3) < ke and \\tp\\ C 2, a(s3) < ne, (35) 

i/ie function 

u(e,T,ip,ip; •) := u x + H l (p,ip; ■ / R e ) + v(e,r, ip, ip ; •), 

solves $25\) in Br c . In addition 

\\v(e,T,(p,tl>; OIIcJ-o-cr*) ^ 2c k£ 2 - (36) 



Observe that the function v{£,T,tp,t()] •) being obtained as a fixed point for contraction map- 
ping, it depends continuously on the parameter r. 
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7 The second nonlinear Dirichlet problem 

For all e G (0, r 2 ), we set 

r £ = y/e. 

Recall that G(x, •) denotes the unique solution of 

A 2 G{x,-) = 64ir 2 S x 

in D., with G(x, •) = A G(x, ■) = on dil. In addition, the following decomposition holds 

G{x, y) = -8 log \x-y\ + R(x, y) 

where y i — ► R{x, y) is a smooth function. 

Given x 1 , . . . , x m G fl. The data we will need are the following : 

(i) Points Y := (y 1 , . . . ,y m ) G Q m close enough to X := {x 1 , . . .,x m ). 

(ii) Parameters A := (A 1 , . . . , A m ) G R m close to 0. 

(hi) Boundary data $ := (cp 1 , . . . , cp m ) G (C 4 ' Q (S' 3 )) m and * := 0\ . . . ,^ m ) G (C 2 ' Q (S' 3 )) m 
each of which satisfies ()22|l . 

With all these data, we define 

m m 

u := ^(1 + A^) G(y\ .) + ^ Xro (.- yi) H e (^,^ ; (• - y*)/r e ) (37) 

3=1 i=i 



where \r Q is a cutoff function identically equal to 1 in B ro i 2 and identically equal to outside 
We define p > by 



B r „ 



P 



4 _ 384 e 4 



(1 + e 2 ) 4 ' 
We would like to find a solution of the equation 

A 2 u - p 4 e u = 0, (38) 

which is defined in Q re (Y) and which is a perturbation of u. Writing u = u + v, this amounts to 
solve 

A 2 v = p 4 e Q+i - A 2 u. (39) 

We need to define an auxiliary weighed space : 
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Definition 7.1 Given r e (0, ro/2), k e M, a e (0, 1) and eel, we define the Holder weighted 
space C^,' a {Vtf (X)) as the space of functions w € C fc:Q (i7f (X)) which is endowed with the norm 

m 

Nlc*.«(fi,(x)) : =\\ w \\c»-«(n ro/2 (x))+J2 sup ( r ""ll«'( a!, '+'-)llc*.'»(fi a -B 1 ))- 

J = 1 r£[r,r /2) 



For all a e (0, ro/2) and all Y e Q. m such that \\X - Y\\ < ro/2, we denote by 

4,y : C°' Q (fi„ (F)) — » C°' a (ft* (Y)), 
the extension operator defined by £ ay y{f) = f in ^o- (Y) 

for each j = 1, . . . , m and £ CT ,y(/) — in each B a / 2 {y :j ), where t i — ► x(£) i s a cutoff function 
identically equal to 1 for t > 1 and identically equal to for £ < 1/2. It is easy to check that 
there exists a constant c = c(v) > only depending on v such that 

l|4,r H|| c o,a (ft , (x)) < c |Hb«(n„ (x))- ( 4 °) 

We fix 

and denote by t/^y the right inverse provided by Proposition l4.il Clearly, it is enough to find 
veC^ a (Q*(Y)) solution of 

v = N(e,A,Y,&,i&) v) (41) 

where we have defined 

N(e, A, Y, $, * ; 5) := £,, Y o £ re , r (p 4 e a+D - A 2 u) . 
Given k > (whose value will be fixed later on), we further assume that $ and \& satisfy 

||$||(C*.«(S3))m <«e, and ||*|| (C 2, Q(S 3 ))m < ne. (42) 

Moreover, we assume that the parameters A and the points Y are chosen to satisfy 

|A|<K£, and \\Y-X\\<K^e. (43) 

Then, the following result holds : 
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Lemma 7.1 Given k > 1. There exist e K > 0, c K > and c K > suc/i i/iai, for all s £ (0,e K ) 

||JV'( e ,A,Y ) *,*;0)|| c j.. (fl . (y)) <c K e 3/2 . (44) 

Moreover, 

||#(e, A, Y, $, * ; i> 2 ) - JV(e, A, F, $, * ; Ui)|| C 4,a (n , (y)) < c K e 2 ||5 2 - Si|| C 4,a (n , (y)) (45) 

and 

\\N(e, A, Y, $ 2 , * 2 ; v) - AT(e, A, Y, $ 1: ^ ; v)\\fU, a(a . (Y)) 

" (46) 

^CkE 1 / 2 (||$2-*l||(C*.°(S3))'" + ||*2-*1 1| (C2,« (S 3 ))m ) 

provided v = v, Vl ,v 2 € C^' Q (f7* (F)), $ = $, $ 1( $ 2 € (C 4 < Q (S* 3 ))"\ # = *,#!, * 2 € (C 2 < Q (S" 3 )) m 
satisfy 

Nlci-(fi*(y)) ^ 2c - £3/2 > ll*ll(c«.-(s»))- < "■ 11*11 (c».-(s»))'» < «£> 

and |A| < ke, ||F- Jf|| < Ky/e. 

Proof : The proof of the first estimate follows from the result of Lemma 15.21 together with l|42|) . 
More precisely, we have 

\\P 4 efi Hc°- (n re (Y)) < c « ^ (4 ^ )/2 and || A 2 u|| c o i(f!re (y)) < c K e 3 / 2 . 

The proof of the first estimate follows from 1J38JI and Proposition l4.il 
The proof of the second estimate follows from 

||p 4 (e a +^ - e^llco-^ {Y)) < c K e 2 \\v 2 - «xlb* (n . {Y)) 

and the third estimate follows from 

IIP 4 (e fi2+11 - e ai+1, )|| c o LVn? . e (y)) < c k£ ( 4 "^/ 2 (||$2 - $i|| (C 4, a(s3))m + ||*2 - *i|| (c2 ,« (S 3 ))m ) 

and 

||A 2 (Q 2 - Ul)ll C °.f 4 (n re (Y)) < C « ffl/2 (11*2 - *l||(C4,a (S8)) m + ||*2 - *l||(Ca.-(S3))-) 

(where iij corresponds to u when $ = <!>., and * = *j) together with (|38[l and Proposition 14. II 

□ 

Reducing e K is necessary, we can assume that 

c K e 2 < - 
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for all e £ (0, e K ). Then, (|44[> and i|45|) are enough to show that 

v i — ► JV(e, A, Y, $, * 5) 
is a contraction from 

{«€C^(n*(y)) : \\v\\ c , a{ ^ {Y)) <2c K e^} 

into itself and hence has a unique fixed point v(e, A, Y, $, ^ ; •) in this set. This fixed point is a 
solution of H39|) . Reducing e K if this is necessary, it follows from l|45|l and H4tj|) in Lemma I7!T1 that 

Me, A, Y, $ 2 , * 2 ;■)-«(£, A, Y, *i,*i;-)ll 

(47) 
< 2c K e 1 ' 2 (||$2 - *i||(c*-»(s3))". + ||*2 - *i||( C 2,o (S 3 ))m ) . 

We summarize this in the : 

Proposition 7.1 Given k > 0, there exists e K > and c K > (Wfa/ depending on k) such that 
for all e £ (0, e K ), /or a/Z sei of parameters A, points Y satisfying 

|A|<«e, and ||Y - X|| < /t\/e 

and boundary functions $ and \£> satisfying \2fy) and 

11*11 (c^(s 3 ))™ < «e, a"rf ll*ll(c 2 -°(S3))" < ke. 
the function 

m m 

u(e, A, Y, *,*;•) := £(1 + A'JG^ + £>,*(• -^V,^' ; (• -y j )/r e ) 

3=1 3 = 1 

+ v(e, A, Y, *,*;•), 



solves fffiSj) m J7 re (K). /n addition 



v(e,A,Y^,^;-)\\ c ^ ( ^ ) <2c K e 3 / 2 . (48) 



Observe that the function t> £l A,Y,*,* being obtained as a fixed point for contraction mapping, 
it depends continuously on the parameters A and the points Y. 

8 The nonlinear Cauchy-data matching 

Keeping the notations of the previous sections, we gather the results of the Proposition 16 . II and 
Proposition l7.ll From now on k > 1 is fixed large enough (we will shortly see how) and e € (0, s K ). 
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Assume that X — (x 1 , . . . ,x m ) £ il m is a nondegenerate critical point of the function W 
defined in the introduction. For all j = 1, . . . , m, we define r» > by 

- 4 logri = R{x i ,x i ) + ^G{x\x 3 ). (49) 

We assume that we are given : 

(i) points Y := (y 1 , . . . , y m ) £ n m close to X := (x 1 , . . . , x m ) satisfying g3J. 

(ii) parameters A := (A 1 , . . . , A m ) £ W n satisfying lEEH) . 

(iii) parameters T :— (t 1 , . . . , r m ) £ (0, oo) m satisfying IJ2l?)l (where, for each j = 1, . . . , m, r* is 
replaced by rl ) . 

We set 

Rt :=r J /x/i 



First, we consider some set of boundary data 

$:=(tp\...,<p m )£ (C 4 ' Q (5 3 )) m and * := (V- 1 , . . . , ^ m ) £ (C 2 ' Q (5 3 )) m 

satisfying {TSJl and (|2~§j) . 

Thanks to the result of Proposition lti.il we can find u, n < a solution of 

A 2 u - p 4 e u = 

in each B re (y J ) , which can be decomposed as 

u lnt (e, T, Y, $, tt ; x) := u e ^ (x - yi) + #V> j ; (x - yi)/r e ) 

+ v(e,T j ,(p j ,tl} j ; Ri(x-y j )/r e ) 

inB P ,(v»'). 

Similarly, given some boundary data 

$:=(£ 1 ,...,0 m )e (C 4 < Q (S* 3 )) m and * := ($\. .. ,^ m ) € (C 2 ' a (S 3 )) m 
satisfying IL'l't and (|42|l . we use the result of Proposition 17. II to find u ext a solution of 

A 2 u - p 4 e u = 
in Q Te (Y), which can be decomposed as 

m m 

u ext (e,A,Y,$,*;x) = ]T(1 + A*) G(y>,x) + X>o(* - V j )H*{(f?^ ; (a; - ^)/r e ) 
+ S(e,A,Y;$,#; a:). 
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It remains to determine the parameters and the boundary functions in such a way that the 
function which is equal to Umt in Uj B Te (y 3 ) and which is equal to u ex t in f2 re (Y) is a smooth 
function. This amounts to find the boundary data and the parameters so that, for each j = 
l,...,m 

Uint = Uext, d r u mt = d r u ext , Au int = Au ext , d r Au int = d r Au ext , (50) 

on dB re (y J ) . Assuming we have already done so, this provides for each e small enough a function 
u e G C 4 ' Q ($1) (which is obtained by patching together the function Ui nt and the function u ext ) 
solution of A 2 u — p 4 e u — and elliptic regularity theory implies that this solution is in fact 
smooth. This will complete the proof of our result since, as e tends to 0, the sequence of solutions 
we have obtained satisfies the required properties, namely, away from the points x 3 the sequence 
u e converges to J2j G(x 3 , ■ ). 

Before, we proceed, some remarks are due. First it will be convenient to observe that the 
functions u 6jT j can be expanded as 

ii e , T i (x) = -8 log |x| - 4 log t 3 +0(e) (51) 

near dB Te . Also, the function 

771 

£(l + A')GV,x) 
i=i 
which appears in the expression of u ext can be expanded as 

m 

£(1 + X e ) G(y e , y 3 + x) = -8 (1 + X 3 ) log |x| + E,(Y; y 3 ) + VE^Y; y 3 ) ■ x + 0(e) (52) 

near dB re (y 3 ). Here, we have defined 

^■(Y;-):=W, -)+E G ^-)- 

In (|50|) . all functions are defined on dB Te (y 3 ), nevertheless, it will be convenient to solve, 
instead of (|50|) the following set of equations 

(u int - u ext )(y 3 +r £ -) = 0, (drUint - d r u ext ){y 3 + r e •) = 0, 

(53) 
(Au int - Au ext )(y 3 + r e ■) = 0, (9 r Au irit - d r Au ext )(y 3 + r e •) = 0, 

on S 3 . Here all functions are considered as functions of z £ S* 3 and we have simply used the 
change of variables x — y 3 + r e z to parameterize 9B re (y- 5 ) . 

Since the boundary data satisfy (|18f) and l|22|l. we decompose 

$ = $ + *i + & X * = 8$o + 12 $i + ^ ± 
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and 

<l = $o + $1 + ^ * = *i+* ± 

where the components of <&o, $ o are constant functions on S 3 , the components of $1 , $1 , ^1 belong 
to Ker(A S 3 + 3) = Spanjei, . . . , e^} and where the components of ^ ± , & x , $ ± , 1 J>- L are L 2 (S 3 ) 
orthogonal to the constant function and the functions e\,...,e&. Observe that the components of 
\P over the constant functions or functions in Ker(A S 3 + 3) are determined by the corresponding 
components of $. Moreover, \P has no component over constant functions. 

We first consider the L 2 (S 3 )-orthogonal projection of (|53fl onto the space of functions which 
are orthogonal to the constant function and the functions ei, . . . , e±. This yields the system 

^,-L_^/,-L = Mo 07 ' 5 (e, A, T,y, $,$,*,§) 

d r ff i (^- L , V- 7 ''- 1 ; - d r H«{<p>> ± ,$> 1 - ; •) = m[ 3) (e, A, T, F, $, $, *, §) 

^7,-L _ ^J,-L = M W ( e , A, T, F, $, $, *, *) 

d r A H i (ip j ' ± ,ip^ ± ; •) - SrAF^.f 1 ; •)) = M 3 (j) (e, A, T, F, $, $, *, §) 



(54) 



where the functions AfP are nonlinear functions of the parameters e, A, F, T and the boundary 
data <&, $, \fr and ^>. Moreover, using (|51|1 and Ij52(l and also l|36|) (keeping in mind that /j, <E (1, 2)) 
and (|48|l (keeping in mind that v G (—1,0)), we conclude that, for each j = 1, ...,m and 
k = 0,1,2,3 

||M fc w) ||c*-» lB(s8) <ce (55) 

for some constant c > independent of k (provided e € (0, e K )). 

Thanks to the result of Lemma f5. 31 and l|55|l . this last system can be re- written as 

($- L , $->-, *- L , vF) = M(e, A, T, F, $, $, *, #) 

where 

l|Af||( C 4, Q ( S 3 )) 2 m><(C 2, Q(S 3))2 m < CE 

for some constant c > independent of k (provided e € (0,e K )). Moreover, (|34J) and l|47|l imply 
(reducing e K if necessary) that, the mapping M is a contraction from the ball of radius ne in 
(C 4 ' a (S 3 )) 2m x (C 2 ' a (S 3 )) 2m into itself and as such has a unique fixed point in this set. Observe 
that this fixed point depends continuously on e, A, T, Y and also on $o, $o, $i, $i and ^i- 

We insert this fixed point in (|53(l and now project the corresponding system over the set of 
functions spanned by e\, . . . , e^ and finally over the set of constant functions. 
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(56) 



The first projection yields the system of equations 

$1 = Mi(e,A,T,y,#o,*o > *i,*i,*i) 

#1 = M 2 (e, A, T,Y,$ ,$ , *!,*!,*!) 

*i = M 3 (e,A,T,Y,*o,*o,*i,*i,*i) 

[ ViV^(F;^') - M 4 W (£,A,T,Y,$ ,$o,$i,*i,*i) 

where the functions Mk (and also M± ) are nonlinear functions depending continuously on the 
parameters £, A, T, Y and the components of the boundary data $0, $o, $i, $i and ^i. 
Moreover, 

|M fc | <ce 

for some constant c > independent of k (provided e € (0, e K )). 

Let us comment briefly on how these equations are obtained. These equations simply come 
from 15U|) when expansions 151|) and (|52|l are taken into account, together with the expression of 
H l (<p>,ipi ; •) and H e ((p> ,ij)i ; •) given in Lemma f5. II and Lemma T5.2I and also the estimates (|36|l 
and l|48|l . Observe that the projection of the term x — ► V£^(Y; y J ) ■ x which arises in Q52JI. as 
well as the projection of its partial derivative with respect to r, over the set of constant function is 
equal to 0. Moreover, this term projects identically over the set of functions spanned by ei, . . . , e 4 
as well as its derivative with respect to r. Finally, its Laplacian vanishes identically. 

Recall that wc have define in the introduction the function 



w(y) :=j2 R (y j >y j )+J2 G(y jl ,y j2 ) 

i=i ji¥=J2 

Using the symmetries of the functions G and R, namely the fact that 

G(x, y) = G(y, x) and R(x, y) = R(y, x) 
we get 



Now, we have assumed that the point X = (x 1 , . 
functional W and hence 

VW\ X 

and 



VW\ Y = 2 (VEi(Y; y 1 ), ..., VE m (Y,y m )). 

,x m ) is a nondegenerate critical point of the 



0. 



(M 4 ) m 3Zh^ D(VW)\x(Z) G (K 4 ) m 
is invertible. Therefore, the last equation can be rewritten as 

,/i(Y-X)=M B (e,A,T,Y,$o,*o,*u*i,*i) 
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The projection of (|53|l over the constant function, leads to the system 

' {\og\/e)- l \og{^ I ri) = M 6 (e,A,T,Y,*o,$o,$i,$i,tti,*i) 

$ = M 7 (e,A,r,y,*o.*o,*i,*i,*i,*i) 

$o = M8(e,A,r,y,#o,*o,*i,*i,*i,*i) 

A = M 9 (e,A > T > y,#o,*o,*i,*i,*i,*i) 

where the function Mk satisfy the usual properties. If we define the parameters U := (u , 
where 

u 3 = — log(r J /r^) 



(57) 



,u m ) 



log 1/e 



and 



Z = y/£(Y-X) 

so that the system we have to solve reads 

(17, A, Z, $ , *o, *i, *i, *i) = M(e, U, A, Z, fc , $ , §i, $i, #i). 



(58) 



where as usual, the nonlinear function M depends continuously on the parameters T, A, Z and the 
functions $o, ^0; 3>ii ^l an d is bounded (in the appropriate norm) by a constant (independent of 
e and «) time e, provided e € (0,e K ). Observe that 

U, A e W n , Z e (M 4 )" 1 , $ io£ K m 

*i,*i,*i £ (Ker(A S 3+3)) m . 

In addition, reducing e K if necessary, this nonlinear mapping sends the ball of radius k e (for the 
natural product norm) into itself, provided k is fixed large enough and e £ (0,e K ). Applying 
Schauder's fixed point Theorem in the ball of radius ke in the product space where the entries 
live yields the existence of a solution of (|58|l and this completes the proof of Theorem ll.il 

9 Comments 

Let us comment on how the condition " (x , . . . , x m ) is a nondegenerate critical point of W" enters 
in our analysis since, we confess, that it is somehow very well hidden. 

The condition " (x 1 , . . . , x m ) is a critical point of W" enters in the estimate ()52fl when Y = X 
and A = 0, since, in this case we have 

m 

Y^ G(x e , x 3 +x) = -8 log \x\ + E 3 {X- x j ) + 0(e) 



2! 



while, if (a; 1 , . . . , x m ) were not a critical point of W, then VEj(X; x^) ^ and we would only 
have 



£<V 



= -%\og\x\+E :j (X-x 1 )+0{e 



1/2N 



which would not be enough : roughly speaking this says that the approximate solution we have 
constructed is not close to any solution of the problem. Given the result of Lin and Wei [7], the 
condition on " (a; 1 , . . . , x m ) being a critical point of W" is a natural one. 

The origin of the " nondegeneracy" assumption takes its roots in the result of Lemma lSTTI which 
classifies all the solutions of the linearized equation about the rotationally symmetric solution. 
The existence of elements 4>i, for i = 1, . . . , 4 in the kernel of L has forced us in proposition 13.11 
to work with weights fj, > 1 to obtain the surjectivity of the operator L M . This choice has one 
importance consequence : In Lemma l5.ll we had to restrict our attention to boundary data which 
satisfy the constraints l|18|) and 11' I'll (even though only the second constraint in (|18|) is important 
to understand where the nondegeneracy condition comes from) to obtain bi-harmonic extensions 
in the unit ball which vanish at the origin at least quadratically. A second reading will convince 
the reader that this property was crucial in the estimate of Lemma 16. II Indeed, the main estimate 
in this Lemma arises from the fact that 

\W{^;-/R £ )\<c^e 2 \x\\ 
Without the second hypothesis in l|18|) we would only have 

\H\^^ : -/R 6 )\<c^e 3 ^\x\ 
which would have led in Lemma l6.1l to the estimate 

||iV(e,r,^,V;0)|| c 4, Q(R)4 <c rt e 3 / 2 

But since /i <E (1,2) this implies that, on the boundary dB^ E the function v(e,T,(p,ip; •) is 
bounded by a constant times e( 3- ^)/ 2 and since 

£ (3-|0/2 >y £ 

the function v would be much larger than the functions H l (ip,ip; -/R £ ) on this boundary and 
hence could not be considered as a small perturbation anymore. Given the fact that, in the 
construction of H % and H e we could not prescribe any function, we had to "find" new degrees 
of freedom to compensate the constraints imposed by (|18fl and l|22|) . The introduction of the 
parameters r J and A J enter at this point to overcome the fist condition imposed by (|18fl and also 
the condition imposed by l|22|l . The points y- 7 close to x^ are introduced to compensate the second 
condition imposed by l|18f) and this is precisely were the nondegeneracy of the critical points of 
W comes into play. 

Let us point out that the nondegeneracy condition strictly speaking can be weakened as this 
has been done for example in [5] and [5] in the case of equation ©. The idea being that the 

25 



nondegeneracy is essentially used to solve the last equation in l|56|) by some disguised version of the 
Implicit Function Theorem. But, remembering that the problem we want to solve is a variational 
problem, this last equation can be rephrased essentially as the gradient of a function W e which is 
defined on il k and which converges (in a sense to be made precise) to the function W as e tends to 
W . Some extra work is needed, but in any case, we could have used some variational technics to 
find critical points of this functional. Since nondegeneracy of critical points is a generic condition 
and in order not to make the exposition of this "nonlinear domain decomposition technic" as 
clear as possible, we have chosen not to follow this route. 

10 Further results 

Modifying very little the previous analysis, it is possible to extend the result of Theorem II. II to 
handle more general equations. We will illustrate this on one example. 

As usual, let us assume that fi C I 4 is a regular bounded open subset and let us choose 
z 1 , . . . , z p £ f2 and a 1 , . . . , a p e (0, +00). We would like to extend the result of Theorem II. II to 
the equation 

( p 

A 2 u = p 4 e u -64n 2 J2 alS z< in n 

u = Au = on d£l. 

Namely, we are still looking for solutions which concentrate at some points ar 1 , . . . , x m € f2, as 
the parameter p — ► and, in order to keep the amount of technicalities as low as possible, we 
will assume that the set of concentration points x 3 and the set of singularities z l are disjoint. 
This problem is very much in the spirit of the work of J2j and [5] even though we do no know 
any applications in physics. On the other end solutions of this problem might be of interest to 
understand constant Q-curvature metrics with conical singularities. 

Setting 

v:=u+ 1 -J2^G{z\-) 

we can rephrase the equation satisfied by u as an equation satisfied by v, namely 



(60) 



< 


p 

A 2 v = p i \\\x-z t \ Aa ' e v in fi 




v = Ai; = on dQ 


This equation is a particular case of the more general problem 


r A 2 u = p 4 Ve u in 




[ u — l\u — on <9rL 



(61) 
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where V : £1 — ► [0, +00) is a smooth function. We are still looking for solutions of this last 
equation which concentrate at some points x , ... ,x m , as the parameter p — ► 0. In order to 
keep the technicalities as low as possible, we will assume that the set of concentration points x 3 
and the set of zeros of V are disjoint. 

As in the introduction, we introduce the functional 

m m 

W(x\...,x m ) :=^ R^ ,x j ) + ^G(x j ,x £ ) + 2^logV(x j ). (62) 

3=1 &i 3=1 

It is easy to check that the result of Theorem 11.11 holds when JIJ is replaced by (|61|) and (|62H 
replaces (@J. We briefly describe the main modifications which are needed to prove this modified 
result. 

Only Sections 6,7 and 8 have to be slightly modified. In Section 6, l|25|l has to be replaced by 

A 2 u = 24e u +e 4 g 

where g is a bounded function (in fact bounded in C°' Q (-Bfj e ) by some constant independent of 
e). It is easy to check that the analysis goes through. The presence of the term e 4 g does not 
alter the estimates of Lemma T6. II and in fact, keeping the notations of introduced in the proof of 
Lemma |6. II we have 

\\e* 9 \\ e ^ (BRe) <ce^. 

The result of Proposition 16 . II remains unchanged. Section 7 applies vertabim and Proposition l7.ll 
is unchanged. 

In Section 8, the main modification due is in the definition of Ui nt - Indeed, for each j = 
1, . . . , m we apply the result of the modified version of Section 6 with 

g=^(A 2 \ogV)(y 3 +e-/r 1 ). 
T j 

This induces in each B Te (y J ) a solution of 

A 2 u = p 4 Ve u 

which can be decomposed as 

u int (x) = u ejT i{x-y°) -logV^z) +H l {tp\^ ; ( x - y 3)/r e ) + v{e, t 1 , if 3 , ^ \ R{(x - y 3 )/r e ). 

The remaining of the analysis of Section 8 remains essentially unchanged once the definition of 
Ej is modified into 

£,(¥;■) :=R(y 3 , -) + ^G(y f , ■) + logV(y 3 ). 

We leave the details to the reader. 
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